Parametrization for the Scale Dependent Growth in Modified Gravity 
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We propose a scale dependent analytic approximation to the exact linear growth of density per- 
turbations in Scalar- Tensor (ST) cosmologies. In particular, we show that on large subhorizon 
scales, in the Newtonian gauge, the usual scale independent subhorizon growth equation does not 
describe the growth of perturbations accurately, as a result of scale-dependent relativistic correc- 
tions to the Poisson equation. A comparison with exact linear numerical analysis indicates that our 
approximation is a significant improvement over the standard subhorizon scale independent result 
on large subhorizon scales. A comparison with the corresponding results in the Synchronous gauge 
demonstrates the validity and consistency of our analysis. 
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I. INTRODUCTION 



Cosniological data from a wide range of sources in- 
cluding type la supernovae [IHll, the cosmic microwave 
background Q , baxyon acoustic oscillations 0], clus- 
ter gas fractions 0, Q and gamma ray bursts 0, 
seem to indicate that at least 70% of the energy den- 
sity in the universe is in the form of an exotic, negative- 
pressure component, called dark energy. While the stan- 
dard ACDM framework is the minimal model that suc- 
cessfully accounts for observations j^, there remain nu- 
merous viable alternatives that also pass current exper- 
imental tests. These alternative models can be broadly 
categorized as quintessence or modified gravity 

models [iB, flsl - fsH , and both categories can, with suffi- 
cient tuning, replicate the expansion history of the uni- 
verse in consistence with observations. (See [s^ [33| for 
recent reviews). In order to distinguish between these 
two categories of models, it is therefore important to look 
beyond the expansion rate. 

The growth of structure offers a hope of breaking this 
degeneracy since different growth histories can arise from 
models which have similar expansion histories. To exam- 
ine the growth of structure one examines the evolution 
of the linear matter density contrast S = 6p/p which 
is given in terms of the background density p and the 
perturbation Sp. For scales much smaller than the hori- 
zon, S satisfies a simple equation called the growth equa- 
tion, which is scale-independent. Note that in what fol- 
lows we use the usual definition of 'scale' as either the 
distance Ap in physical FRW coordinates or the corre- 
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sponding wavenumber k = j^. There are many inves- 
tigations which attempt to characterize the evolution of 
6 through the use of a growth parameterization which 
assumes a different value depending on the cosmological 
model used, thereby allowing for models to be distin- 
guished (e.g. [H, d^). The standard definition of 
the growth parameter, 7, in terms of the growth func- 
tion /, the matter density flm and the scale factor a is 
given as 
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Once this parameter is determined (for earlier theoreti- 
cal developments on the parameterization of the growth 
parameter and experimental constraints on 7, see (33.1361- 
|47[) one may then be in a position to determine whether 
the standard general relativistic (GR) framework of 
ACDM is responsible for the acceleration of the universe, 
or some other, more exotic process is at work. 

However, in [H, [i^, it was demonstrated that in the 
Newtonian gauge (for another interesting look at gauge 
issues see [50[) the above parameterization can become 
inaccurate for large subhorizon scales > 100/i~^Mpc. In 
other words, if the physical growth of structure is cor- 
rectly described by the Newtonian gauge, then it would 
show up as inconsistent with scale-independent param- 
eterization in Eq. ([T]) and (mistakenly) appear to be 
caused by exotic physics. The reason for the discrepancy 
was shown to be the scale dependence of the growth of 
6, which becomes important for large subhorizon scales 
(> 100ft.~^Mpc). An improved version of the growth 
equation was derived in |48| which incorporates the scale- 
dependence. In [49j . a new scale-dependent parametriza- 
tion of the growth function / was proposed, which was 
shown to account for the evolution of /(a) on these large 
scales with considerably greater accuracy than Eq. ([Ij. 

In the present work, we focus on the growth of per- 
turbations in scalar-tensor (ST) theories of gravity. It 



2 



is well known that in the sub-Hubble approximation the 
growth of perturbations in these theories is also described 
by an equation similar to the growth equation in GR, up 
to a redefinition of the gravitational constant. Working 
in the Newtonian gauge, we show that the usual growth 
equation approximation for S becomes unreliable on large 
scales, for the same reason as in the GR case, i.e. the ef- 
fects of scale dependence. We derive an improved version 
of the growth equation relevant for these models and pro- 
pose a more accurate scale-dependent parameterization 
for growth. 

The layout of our paper is as follows. In Section 
II we discuss the growth of perturbations in ST theo- 
ries, demonstrating the failure of the usual growth equa- 
tion approximation and introducing an improved growth 
equation and a new parameterization for the growth func- 
tion in these models. In Section III we compare these 
approximations to exact solutions to demonstrate their 
accuracy. Our conclusions can be found in Section IV. 



II. GROWTH OF MATTER PERTURBATIONS 
IN SCALAR-TENSOR GRAVITY BEYOND 
SUBHORIZON SCALES 

ST theories are widely studied as an alternative to 
GR. These theories arc well-motivated from string the- 
ory, Randall-Sundrum models, as well as extended and 
hyperextended inflationary models. (See e.g. [5l| and 
references therein for a review of ST theories.) The de- 
viations from GR predicted by these theories have been 
investigated (see e.g. j53 - l59| ). ST theories have also been 
used to explain the accelerating Universe, and the cos- 
mological consequences of these models have been widely 
studied (see e.g. [13, [H, UMl)- 

In this work, we focus on the growth of matter per- 
turbations in these theories. As shown in [13,13' if o^^^ 
works in the Newtonian gauge and considers scales much 
smaller than the horizon (k 3> aH) then the overdensity 
6 obeys an equation very similar to the familiar 

6 + 2H6-'lTrG,s{t)p,„S = 0, (2) 

where dots denote derivatives with respect to cosmic time 
and Gesit) is an effective gravitational constant whose 
evolution is determined by the scalar field dynamics (see 
also equation below). We reconsider the growth of 
perturbations in these models, manifestly retaining the 
scale-dependent effects, and derive an improved version 
of the growth equation which models the evolution of S 
with a greater accuracy than the scale-independent equa- 
tion. Using our improved growth equation, we then pro- 
pose a new scale-dependent parameterization for growth 
in these models which is applicable under the assump- 
tion of our approximations which involve slow evolution 
of the scalar field. 



A. Scalar-tensor cosmology 

We start with the general action of a Universe de- 
scribed by ST gravity (in the Jordan frame) and arbitrary 
matter fields: 

S - ji^J d^^V-9[Fim (3) 

where y^j^ is the metric with determinant g and Ricci 
scalar R. G* is the bare gravitational coupling constant 
(henceforth we will set SttG* = 1). The scalar field 3> has 
a potential U (<&) and couples to gravity through the func- 
tions F(^) and Sm denotes the action of matter 
fields ijjfn ■ Henceforth we work with the paramctrization 
where Z($) = 1 and -F($) is arbitrary. 

We next consider a spatially fiat Friedman-Robcrtson- 
Walkcr (FRW) Universe with a background metric (in 
the Jordan frame): 

ds^ = -dt^ + a^{t) [dr^ -f [dO^ + sin^ OdLp"^)] . (4) 

We take the matter content of the Universe to be a perfect 
fluid with energy-momentum tensor: 

Tf,^ = -^^^ = {p + p)Uf,u„ +pgf,„ , (5) 

where p, p and arc the pressure, energy density and 
four-velocity of the matter fluid respectively. For con- 
venience, we will henceforth consider the matter to be 
pressureless and set p = 0. 

Finally, we work at linear order in perturbation the- 
ory in the Newtonian gauge. The metric is perturbed as 
follows: 

ds^ = - [l + 24)) dt^ + ci^ [l - 2^:) dx^ (6) 

where (j) and ip are the temporal and spatial perturbations 
and X = {r,9,ip). The matter energy density is perturbed 
as p — > p + where p is the background density and 
6p is the perturbation. For convenience, we define the 
overdensity S = Sp/p. The perturbations in the velocity 
field Sufi are conveniently expressed through the velocity 
potential v (defined such that (Su^ = —d^v). The scalar 
field is perturbed as $ — > <!> S^. 

The evolution of the background variables is governed 
by the zcro'th order Einstein equations and the equations 
of conservation of energy-momentum: 

3FH^ = + - iHF + U (7) 

- 2FH = p.m + ^^ + F - HF (8) 
<i> + 3H^ = 3F'{H + 2H^)-U' (9) 
p„, + iHpra = (10) 

where primes denote derivatives with respect to the scalar 
field $ and dots denote derivatives with respect to the 
coordinate time. 
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The evolution of the perturbations is governed by the 
first order Einstein and conservation equations. We 
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work in Fourier space assuming a spatial dependence of 
exp(ik.x). The scalar field fluctuation is given by 



(5$ + 3i7(5$ 



-3{H + 2H^ F" + U 



(5$ = 



2V') - 3 (i/i + mi) + 



F' 



(11) 



r 



The matter density perturbation and velocity potential 
evolve as 



i) = (j) 



(12) 
(13) 



The evolution of the metric perturbations is given by 
the equations 



-3F<j>- 



ip = (/) + 



2F 



2^F- $2 



(^ijj + H(i^ + Fc/) 
3HF 



F'S'^ 



F 



pv + ^S^ + SF - HF'S^ 



p{5 + iHv) + C/'(5$ 



F ■ 

+3iJi> 5^ + 3-5F 



F2 



2 QH'-rj^]SF + ^S'S> 



(14) 
(15) 



(16) 



r 



B. Growth of matter perturbations 



Note that Eans. (|12ll3p can be combined to obtain an 
exact second-order differential equation for S: 



• fc2 

6 + 2HS +—(l)^3Th + 6H1P 



(17) 



We now proceed to obtain an approximate equation for 
the growth of perturbations of subhorizon modes. How- 
ever, instead of considering only terms which are pro- 
portional to fc^/a^, as is usually done (see e.g. [isl. Is^). 
we also retain terms proportional to H^. This makes 
the analysis more accurate for sub-Hubble modes close 
to the Horizon size. In addition, we also assume 
matter-domination, during which the metric potentials 
are frozen. Thus we ignore time derivatives of metric 
perturbations. We also ignore time derivatives of the 
field. Even though this approximation is not always ap- 
plicable, in practice we have found that our analytical 
scale dependent approximation for the growth of pertur- 
bations turns out to always be significantly better than 
the usual scale independent sub-Hubblc approximation. 



The precise accuracy level however depends on the de- 
gree of validity of the above approximation which may 
vary depending on the details of the field dynamics. 
Defining 



3-ffo^o 



(18) 



(where the last equation is valid in the matter era) 
Eq. pT|) reduces to: 



(5$ 



FF' + 2^FU'/3H^ 



F + 2F'2 + ^ {FU"/3H^ - F"U/2H^ - FF"/2) 

(19) 



Plugging into Eq. (fT6|) we obtain 



'^(b=-^piS + 3Hv)g{F,U,0 , (20) 



where 
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2F + 4F'2 + ^ {2FU"/3H'^ - FF" - F"U/H^ 



2F + 3F'2 + ^ (2F'2 - FF" + 2FU" /3H'^ - U'F'/H^ - F"U/H^) + (AF'U' /3H^ - 2U'^/9H^) 



,p6 



2^^ \l + ^9iF,U,0h{F,U,0 



(21) 



r 



Using the above, (and eliminating v from Eq. psp ignor- 
ing time derivatives), we obtain the foUowing form for 
the Poisson equation: 



where 



h{F,U,0^ 



2F + 3F'2 + ^ (2FU"/3H^ - 2F'U'/3H^ - F"U/H^ - FF") 
2F + 4F'2 + ^ (2FC/"/3iJ2 - F"U/H^ - FF") 

I 



This leads to our "improved" growth equation for the 
evolution of perturbations, aecurate for large subhorizon 
scales: 



6 + 2H6-± 



1 + ^9 iF,U,OhiF,U,0 



pS = 0. (22) 



In the case of general relativity this reduces to the form 
derived in Ref. HEl 



5 + 2HS 



(23) 



Note that for scales much smaller than the horizon 
(fc 3> aH, or equivalently ^ — > 0) equation 
to the well-known form (see e.g. Eqn.(5.13) of 



2Hd - — 
2F 



2F + 4:F'^ 
2F + 3F'2 



pS = 0. 



(24) 



Equation ()23p may be expressed in terms of the growth 
factor f — iu the form 

*' din a 



f' + f+(2-Ma)]f=- 



3 fim(a) 



2 1 + e(a,fc) 



where ' = ,f , and we have assumed ACDM for H(a). 

For sub-Hubble scales £,{k,a) — > and the solution of 
equation ((25)) is well approximated by ([1]) with 7 = ^■ 
For larger scales a perturbative approach [i^ may be 
used to derive the scale dependent growth rate /(a, k) of 
matter perturbations in GR as: 



(26) 



with K = for the scale-independent growth function 
and K = 1 for the GR scale-dependent approximation. 



Following the same reasoning as in Ref. [431, the fol- 
lowing parametrization may be derived for the growth 
hmction f{a,k) in ST cosmologies: 



f{a,k)^n,niar- 



9{F,U,0 



l + ^g{F,U,Oh{F,U,C) 



(27) 



Even though this parametrization may be derived as an 
approximate solution by a perturbative expansion to all 
orders in ^ as in Ref. [4^, its validity for 7 ~ 0.55 (the 
exact value is obtained by fitting to the exact numerical 
solution) will also be tested in the next section. The 
scale dependence of 7 has also been recently discussed in 



III. COMPARISON TO EXACT RESULTS 

We now compare the evolution predicted by Eq. ([^^ 
and Eq. p7)) to results from the exact evolution. 

To solve the system numerically, we choose the follow- 
ing functional forms [85|. We take 



(25) ^($) = 1 _ Xf^^ ; Z($) = 1 ; L/($) = 1 + exp[-A$] 

(28) 

Then we solve numerically the background equations 
given by Eq. ©-Eq. (HUl) (for details see [ssj). 

To evolve the perturbations, we solve equations 
Eq. (US]), Eq. (HH) and Eq. ([TTI). These are solved in 
the Newtonian gauge using the background. The initial 
time is taken to be in the matter dominated epoch, at 
Zi = 1000, and we set $i ~ 0. 

In our numeric solution the initial value is cho- 
sen so that at early times the deviation in the expan- 
sion rate from the ACDM one is not larger than around 
5%. This deviation is roughly given by ^st-Hacum ^ 



F($i) ^1"^ — which then sets an upper bound 'I'i.max for 
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FIG. 1: For three choices of scale, plot of the exact 5m{k,a) 
(blue solid line), along with two approximations: (i) the so- 
lution of the scale-dependent ST growth equation (|22[) (solid 
line) and (ii) the solution of the scale-independent ST sub- 
Hubble growth equation (dot-dashed line). The apparent 
scale dependence of the present time value of 5 in the sub- 
Hubble approximation is due to the scale dependent initial 
conditions considered needed to secure a scale independent 
value of the initial metric perturbations [sHl 



$. Although $i/$i,max ~ 1 results in an expansion rate 
consistent with the ACDM one, we find that such a choice 
gives rise to deviations between the analytic approach in 
Eq. (22) and the numeric solution for Sm- This is because 
larger $ implies larger effective potential energy of the 
field, and hence larger $ when the field oscillates due to 
its coupling to curvature and therefore to matter. Above 
a certain threshold, $ grows large enough so that one 



FIG. 2: Plots of the % difference between the exact numer- 
ical solution for Sm (blue line) and the scale-dependent ST 
approximation (solid line). The % difference corresponding 
to the subhorizon ST approximation is not included since this 
is much bigger than the corresponding to the scale-dependent 
approximations. The values used for the plot are A/ = 5 and 
A = 10. 



cannot keep neglecting it in order to arrive to Eq. 
In what follows we consider $i/$i,max ^ 0.30. For this 
initial condition we secure relatively small deviation from 
GR and from the ACDM expansion rate (consistent with 
observations) and small field time derivative (consistent 
with our approximation) . We thus find that the solution 
to Eq. (p2)) deviates from the numeric solution for by 
less than 1% on subhorizon scales. 

We first compare the exact numerical solution of the 
ST perturbation system ((TS|) . (potentials given by 
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FIG. 3: Plots of the % difTerence between numerical solution 
for S (in blue) and (i) the scale-dependent ST approximation 
(solid line) and (ii) the subhorizon ST approximation (dot 
dashed line) as a function of the perturbation scale Ap = ^ . 
We set A/ = 5 in both plots and A = 10 and A = 1 for the 
top and bottom panels, respectively. 



value of i5. Clearly the scale dependent approximation 
is in good agreement with the exact solution on physi- 
cal scales up to Ap = 10^/i~^Mpc (Ap = ^). On the 
other hand, the sub-Hubble scale independent approxi- 
mation (dot-dashed line in Figure [Sj starts deviating sig- 
nificantly from the exact solution already on scales larger 
than about 10^/i~^Mpc. On scales significantly larger 
than 10'^/i~^Mpc, both approximations break down as 
shown in Figure |3l 

The approximate scale dependent growth rate f{a,k) 
obtained in equation (|27p is compared with the corre- 
sponding exact result in Figure S] where we show a su- 
perposition of the growth factor /(a, k) corresponding to 
the approximations (sub-Hubble scale independent ST 
and scale dependent ST) and to the exact solution. The 
parameter values used in Figure 2] are A = 10, Xf ~ 1 
which provide a background evolution similar to the best 
fit KCDM model (the difference for the expansion rate 
is always less than 5%). Clearly, the scale independent 
growth factor deviates significantly from the exact result 
while the scale dependent approximation remains very 
close to it even on large scales [k = O.OOl/i Mpc~^) and 
rcdshifts. Thus even though on relatively small scales 
(fc > 0.01/iMpc~i) or low rcdshifts {z < 2) the accuracy 
of the two approximations is good and similar, on larger 
scales and rcdshifts the accuracy of the scale dependent 
approximation is significantly better than the scale in- 
dependent approximation which fails to approximate the 
exact numerical solution. 



IV. CONNECTION TO THE SYNCHRONOUS 
GAUGE 



p8| ) with the solution of the approximate linear growth 
equation assuming the following approximations: 

• Scale Dependent ST growth equation ([22]) . 

• Scale independent ST sub-Hubble approximation 
(HI 

This comparison is shown in detail in Figures [T] - 131 In 
Figure 1 we show the evolution of the matter density 
perturbation S{z) in the above three cases (two approxi- 
mations and the exact solution). For intermediate scales 
{k = 0.01/i~^Mpc, upper panel of Figure [1]) the scale 
dependent approximation of equation p2|) (black con- 
tinuous line for ST) fit well the exact numerical result 
(blue continuous line) while the scale independent sub- 
Hubble approximation (dot-dashed line) is significantly 
less accurate. On larger scales, the sub-Hubble approx- 
imation becomes even less accurate (lower panel of Fig- 
ure [1} while the scale dependent ST approximation of 
eq. (|22p remains a good approximation to the exact nu- 
merical solution within a few percent (Fig. [2]). This is 
more clearly demonstrated in Figure [3] where we show 
the scale dependence of the deviation of the two approx- 
imate present values of S from the corresponding exact 



The above calculations have been made in the confor- 
mal Newtonian gauge. In this section we briefly discuss 
the gauge dependence of the density perturbation S. Let 
us consider the gauge invariant quantity 

5™ = + 3Hv (29) 

P + P 

where v is the velocity potential. This quantity may 
be evaluated in the synchronous and in the Newtonian 
gauges leading to 

Sis) ^ (J^)(^) . (^)W + (3i7.)W (30) 

P+P P+P 

where the superscripts and imply evaluation in the 
Synchronous and in the Newtonian gauge respectively. 
Assuming matter domination {p = 0) and dropping the 
index for the Newtonian gauge we find 

5''^^=6 + 3Hv (31) 

The evolution equation for matter perturbations in the 
synchronous gauge is 

+ 2i7 + 47rG(5'^) = (32) 
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FIG. 4: For three choices of scale, plot of the exact growth 
function f{k,a) (green, solid line), along with two approxi- 
mations: from the scale-dependent ST growth equation (|27|) 
(black, dashed line) and from the scale-independent GR sub- 
Hubble growth equation (|26p with K = (red, dot-dashed 
line). 



which is scale independent. The corresponding equation 
in the Newtonian gauge is eq. (|17p. When properly re- 
lated initial conditions are used in the solutions of eqs. 
^T7\ and ([5^ . their solutions should satisfy cq. (f5T]) . 

Consider for example initial conditions of the form 



-10 



-5 



const 



(33) 



It is straightforward to show (see Appendix) that the 
above initial conditions in the matter era transform to the 
following initial conditions for the matter perturbations 
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FIG. 5: For a large scale {k = O.OOlhMpc'^) we show the 
evolution (with parameter values A/ = 5, A = 10) of S{a) 
(blue continuous line), S^^\a) (black dot-dashed line) and 
S{a) + 3Hv (red dashed line which coincides with the line of 
5(^'(a) as anticipated). 



in each gauge: 



AS) 



2P 



for the synchronous gauge and 

2fc2 



6^ = 

s,. = 



_2P_ 
3^ 



+ 2 V' 



(34) 
(35) 

(36) 
(37) 



for the Newtonian gauge. Using these initial conditions 
we solve the perturbation equation in each gauge (eqs. 
(IT71) and ^) and verify the validity of eq. This 
is an additional verification of the validity of our analy- 
sis and demonstrates that the difference of SHv between 
the matter density perturbations in the two gauges can 
be significant. It is demonstrated in Fig. [5] where we 
consider a large scale (fc = 0.001 h Mpc'-^) and show 
the evolution of (5(a) (Newtonian gauge, blue continuous 
line), (5'^'^'(a) (Synchronous gauge, black dot-dashed line) 
and (5(a) -\-3Hv (red dashed line which coincides with the 
line of S^^\a) as anticipated). 



V. CONCLUSIONS 

We have derived and tested numerically a simplified 
ordinary differential equation whose solution describes 
fairly accurately the growth of linear cosmological pertur- 
bations in ST cosmologies up to Hubble scales (beyond 
10'^/i~^Mpc). A corresponding analytic form of the scale 
dependent growth rate /(a, k) was also presented and 
tested numerically in comparison with the exact linear 
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result. This is a significant improvement over the pre- 
viously known sub-Hubble which breaks down on scales 
larger than about 10^/i~^Mpc. 

We have thus demonstrated that in the context of the 
Newtonian gauge, the comparison of cosmological large 
scale structure data with corresponding theoretical pre- 
dictions in ST based cosmologies on scales larger than a 
few hundred Mpc should not be based on the sub-Hubble 
scale independent approximation. Instead it requires ei- 
ther use of the full linear numerical solution of the cos- 
mological perturbation equations in each theory or the 
imposed scale dependent approximation presented in the 
present study. Thus, previous work using the scale inde- 
pendent sub-Hubble approximation to make predictions 
for the growth of perturbations in ST theories is reliable 
on scales up about 200h''^Mpc. However, such results 
should be used with care on larger scales. For example, 
on scales of a few hundred Mpc the error induced in the 
growth rate by not using the scale dependent effects is 
about 10% on redshifts z > 3 while on scales larger than 
IGpc the corresponding error exceeds 50% even at low 
redshifts z ~ 1. 

We have also shown that our results are consistent with 
a corresponding calculation in the synchronous gauge and 
verified that the quantity ^ is gauge dependent and 
there can be a significant difference between its forms 
in different gauges on large scales. 
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Appendix A: Growth in synchronous and conformal 
Newtonian gauges 

In this appendix we discuss the appropriate (scale- 
dependent) initial conditions that must be used in the 
Newtonian and synchronous gauges in order to verify 
Eq. ([31]), i.e., equations (p4| - (|37| . 

Starting with Eg. p^ . and using the initial conditions 



S^, = (5$, = i', = 0, ~ -10-5 ^ (Al) 

(a subscript i will always refer to the value of a quantity 
at the initial time) we obtain the following relation for Si 




Using Ea. ([T5t with the above initial conditions (|A1[) 
we arrive at the relation 



Next combining Eqns. (|A2[) and (|A3[) . and ignoring 
time derivatives, we can then recover the initial condition 
for S: 

^,--(2^ + 6i.f)^ (A4) 

\ a J pmi 

From Eq. (|X3)) . using Eqs.([ni), jTS]), and ([H]), as weh 
as the initial conditions (|A1[) . we find 

5,^-%— (2H,F, + (A5) 

Also from Eqs.Q and ([8]), ignoring time derivatives, 
we have the approximate relationship 

^^p,^., ^ -^Hf (A6) 

Lastly, ignoring time derivatives in Eqns. (|A3[) and 
(|A5|) . and using Eq. (|A6p . we obtain equations (p4| and 
(ESI). 

For the synchronous gauge, we start with (Eq. ([31]) ). 
use equations (|A2[) and (|A6p , and ignore time derivatives 
to obtain Eq. ([55]) . Then from Eqn. (IT^ . using equations 
(|A1[) . (jA3p . (|A5[) and (jA6p . and ignoring time derivatives, 
we obtain Eq. ([57]) . 
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